Structural stability of the fullerene against defect formation 
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Using a systematic search algorithm we identify several types of point defects in the boron fullerene 
with 80 atoms. All these point defect leave the cage structure intact. In addition the cage structure 
is also very stable with respect to elastic deformations and addition or removal of atoms. 



Since boron is a neighbor of carbon in the periodic ta- 
ble, they share a lot of similarities. The sp'^ hybridization 
of the valence electrons, the electron deficiency, large co- 
ordination numbers and short covalent radius allow boron 
to form strong directed chemical bonds and give rise to 
rich variety of possible nanostructures. Boron nanos- 
tructures share some common basic units following the 
so-called " Aufbau principle" [l| . Various kinds of boron 
nanostructures, such as boron nanotubesjl, boron 
fullerenes [1, [5[ and boron sheets [2] have been studied 
both experimentally and theoretically. The B^o fullerene 
has recently been identified in density functional calcula- 
tions as an energetically very favorable structure [6]. Like 
the carbon fullerenes boron based fullerenes might thus 
play an important roleQ if one succeeds in synthesiz- 
ing them. In contrast to other fullerenes where defects 
have been studied in detail, the defects in have not 
yet been studied. Since materials are never perfect in 
reality, it is important to have information about the 
structure and energetics of their defects, in this article 
we study the defects of Bso- Whereas in previous stud- 
ies of defects in various materials educated guesses were 
made to identify defects, we use a systematic search al- 
gorithm [10] to find them. As a consequence we find 
not only one type of defect but several of them. All the 
defects we found leave the cage structure intact. This 
stability should make their synthesis easier. If in an ini- 
tial stage of the synthesis an imperfect structure is gen- 
erated, it forms nevertheless a cage structure and it can 
later anneal into the perfect cage structure. 

In the first publication on the ^80 fullerene a per- 
fect icosahedral structure was proposed Q. In a sub- 
sequent publication a distorted icosahedron was pre- 
sented as the ground state if the PBEfsj density func- 
tional is used. Our calculations confirm that the dis- 
torted cage is the ground state within PBE. However the 
distortion is marginal and atoms are displaced by less 
than 0.035 Bohr with respect to the ideal structure as 
shown in Fig[Tl The energy differences are also small and 
as a matter of fact a LDA functional [9] favors the per- 
fect icosahedral structure by about the same amount of 
energy (Table [Tj). Because the deviations from the icosa- 
hedral structure are so small we will call in the following 
both the perfect and the slightly distorted icosahedral 



states just icosahedral. 




FIG. 1: A section of the ^so fullerene that illustrates 
the differences in the ground state geometries of the LDA 
ground state (grey spheres) and the PBE ground state (black 
spheres). In the LDA geometry the hexagon centered atoms 
are exactly in the plane of the hexagon, in the PBE geometry 
they move slightly inside. 



In order to search in a systematic way for defect struc- 
tures we used the minima hopping algorithm to- 
gether with the BigDFT electronic structure program 
from the ABINIT package [12]. The BigDET electronic 
structure program uses a wavelet basis set. This is a sys- 
tematic basis set and the number of basis functions was 
chosen such that energies are accurate to within better 
than 1 mHa. Pseudopotentials were used to eliminate 
the core electrons [l3|. Minima hopping is an algorithm 
to find the global minimum as well as other low energy 
structures. It is entirely unbiased and applicable to any 
system. After exploring some 20 configurations with the 
minima hopping algorithm we could not find any con- 
figurations where the cage structure was destroyed. The 
configurations we found in this way were combinations of 
several of the basic defect structures that we will discuss 
below. Combining these basic defects can give rise to a 
large number of defective cage structures all of which are 
very close in energy to the ground state. 

A ^80 fullerene can be obtained from the Ceo fullerene 
structure by adding 20 atoms in the center of the 
20 hexagons of the Ceo fullerene structure. The ^80 
fullerene exhibits however a larger variety of defect struc- 
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tures than the Ceo fuherene. There are defects where a 
pentagon is destroyed and a distorted hexagon or hep- 
tagon is created. We will in the following call this new 
class of defects hexagon/heptagon defects. Once a defect 
has been created it can be shifted to different places on 
the surface of the cage. This is shown in Fig. [2l The 
barriers for the displacements of the hexagons are low. 
For the process leading from the perfect structure (top 
structure in Fig. [2|) to the first defect structure (below 
perfect structure in Fig. [2]) the barrier height is .073 Ha 
(2.0 eV) in LDA and .067 Ha (1.8 eV) in PBE. Displace- 
ments of defects are thus expected to give rise to some 
plastic deformation behavior of the ^80 fullerene. 

Hexagon defects can also be created by shifting a group 
of atoms as shown in Fig. [3l In contrast to the previous 
hexagon defect shown in panel D of Fig. [21 where the 
hexagon is neighbor to a pentagon, the hexagon is iso- 
lated from the pentagons in this case. The fullerene cage 
is also relatively strongly distorted for this defect, we 
will call this defect type 4 defect. 

It is not too surprising that the B^o fullerene inherits 
also the Stone- Wales type of defect |1J|. If one adds to 
a Ceo with a Stone- Wales defect again 20 atoms into the 
hexagons one obtains the ^80 with a Stone- Wales defect 
shown in Fig. [H The number of pentagons is 12 as in 
the perfect structure, but two pentagons became now 
neighbors. More pathways exist however in ^80 than in 
Ceo to create a Stone- Wales defect. One can for instance 
transform the defect in Fig. [3] into a Stone- Wales defect 
by moving one atom which is part of the hexagon and 
closest to the next pentagon towards the center of the 
hexagon. 





LDA 


PBE 


perfect Icosahedron 


0. 


0. 


distorted Icosahedron 


.020 


-.020 


defect type 1 


.053 


.027 


defect type 2 


.062 


.035 


defect type 3 


.027 


.066 


defect type 4 


.061 


.040 


Stone Wales type defect 


.050 


.026 



TABLE I: The energies of various configurations with respect 
to the highest symmetry icosahedral fullerene as calculated 
with the LDA [3] and PBE [8] density functionals in unit of 
Ha. All the energies are the energies of relaxed structures 
with the exception of perfect icosahedron in PBE and the 
distorted icosahedron in LDA. 

In addition to being stable against defect formation, 
the ^80 fullerene is also highly stable against elastic de- 
formations. We enforced for instance several strong de- 
formations onto the cage such as a sharp edge together 
with flat regions shown in Fig. [5l Such a structure could 
serve as a nucleation site for a more compact cluster 
structure. We found however that the fullerene always 



relaxed back to the open icosahedral structure during a 
geometry optimization. 

The ^80 fullerene is also stable against the removal 
or addition of a single atom. If one removes a boron 
atom from the center of a hexagon to obtain Bjq or if 
one adds a boron atom in the center of a pentagon to 
obtain ^81 one obtains structures whose cohesive energy 
per atom differs by less than .5 mHa from the values of 
the icosahedral Bso- Consecutive hexagons can be filled 
until one obtains the Bg2 fullerene. The cohesive energies 
per atom changes hardly when going from ^80 to Bg2- 
The cohesive energy for Bg2 that we calculated is only 1 
mHa less than the one for B^o- Our number is exactly in 
between the two numbers given in reference [6]. 

In summary, we have shown that several different types 
of point defects exist in the ^80 fullerene and that the 
cage like structure is stable against all these point defects 
as well as against strong deformations and changes in the 
number of atoms. Since other boron nanostructures are 
built up according to the same basic construction prin- 
ciples, it is to be expected that the same or very similar 
defects can be found in other boron nanostructures. 
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FIG. 2: Illustration of the hexagon and heptagon defects. 
To get from the icosahedral ground state configuration (top 
panel A) to the hexagon defect structure (panel below) one 
has to shift the black atom in the uppermost panel upward 
as indicated by the arrow. The hexagon defect in the second 
lowest panel B can then be shifted to the positions in the 
second lowest and lowest panel C and D by moving upwards 
or upwards to the left either of the two black atoms in the 
second uppermost panel as indicated by the two arrows. In 
panel B and D the defect is a hexagon, in panel C a heptagon. 
The defects in panels B,C and D will be denoted in Table U 
as type 1,2 and 3. 



FIG. 3: The picture shows the most distorted part of a type 4 
defect in ^so- The hghter atoms and dotted bonds show the 
positions of atoms in perfect structure and the dark atoms 
with sohd bonds shows the Bso with the defect. 




FIG. 4: A perfect Bso (on the left) and a Stone- Wales type 
defect in Bgo- 




FIG. 5: A Bso fullerene that was deformed such that it has 
a sharp edge and two flat surfaces. This strongly deformed 
structure relaxes back to the icosahedral structure. 



